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We show how to detect entangled, bound entangled, and separable bipartite quan- 
^^ ' turn states of arbitrary dimension and mixedness using geometric entanglement wit- 

^ ■ nesses. These witnesses are constructed using properties of the Hilbert-Schmidt 

geometry and can be shifted along parameterized lines. The involved conditions are 
D , simplified using Bloch decompositions of operators and states. As an example we de- 

^^ \ termine the three different types of states for a family of two-qutrit states that is part 

Q> ' of the "magic simplex", i.e. the set of Bell-state mixtures of arbitrary dimension. 



INTRODUCTION 



-)— > 

2 \ Entanglement is a fascinating curiosity of quantum physics that distinguishes it consid- 

Q-i| erably from classical concepts [Ij. On the one hand it implicates surprising philosophical 

aspects such as the incompatibility of local realistic theories with quantum physics [2|, y], on 

CN ■ the other hand it can be successfully implied in quantum information and quantum commu- 

^ ! nication tasks to improve quantum protocols with respect to classical ones (for an overview 

(^ ; see, e.g., Refs y, a, l6|). 

QO ' It is still an open mathematical problem to determine whether a quantum state is en- 

tangled or not, there is no operational procedure for an arbitrary state, although there 
are many useful criteria that allow a detection of entanglement in many cases [71, |8|, |9(]. 

00 ! For pure states and lower dimensional bipartite systems (e.g., two qubits) the problem is 

O ■ solved, since there exist applicable necessary and sufficient conditions for separability (i.e. 

non-entanglement) [10]. Additionally, entangled states can be classified according to their 

KA ', distillability: A distillable state can be "distilled" to a (nearly) maximally entangled state 

H I via statistical local operations and classical communication (SLOCC). States that are not 

distillable are called bound entangled, whereas distillable states are called free entangled 
[Sl, lul, bjl, [l3|. Examples of bound entangled states and construction procedures can be 



found in Refs. 14, 15,hj,hj, 18 



In this article we want to present a method to classify entanglement with entanglement 
witnesses, which provide a well-established tool to detect entanglement. Our approach is 
based on two concepts: We construct witnesses in a geometrically intuitive way and use Bloch 
decompositions of operators and states to simplify the mathematical application. We show 
how to use geometric entanglement witnesses to detect more entangled states than the PPT 
criterion and how to identify separable states by using optimal entanglement witnesses. This 
is attained by shifting the witnesses along parameterized lines of states. Two main methods 
are explained in detail: the outside-in shift and the inside-out shift. The outside-in shift 



Electronic address: philipp. kranimer[at]univie.ac. at 



is used for detecting more (bound) entangled states, whereas with the inside-out shift we 
construct the shape of the set of separable states. 

The paper is organized as follows: In Sec. [TTl we give an overview of the mathematical 
basics of entanglement theory that will be used throughout this article. We present a 
formulation of the entanglement witness criterion in terms of Bloch decompositions in Sec. lIIII 
and explain a geometric method to construct entanglement witnesses in Sec. IIVI where we 
also introduce the two shift methods. The results are applied to a family of states that 
are part of a special simplex (called "magic simplex" ) in the Hilbert- Schmidt space of two 
qutrits in Sec. |Vl where illustrative geometric pictures are obtained identifying regions of 
entangled, bound entangled, and separable states. 



II. BASIC CONCEPTS OF ENTANGLEMENT THEORY 

We consider a bipartite Hilbert-Schmidt space A := Aa ® Ab on a discrete finite dimen- 
sional Hilbert space TC := Ha ® Hb of dimension dA x ds, D := dAdB- Vector states \iIj) 
are elements of 7-^ whereas density operators (called "density matrices" or just "states") are 
elements of A^ a Hilbert space of operators on 7i, with a scalar product 

(A, B) = TiA^B, with A, B E A. (1) 

States p E A are defined by the properties 

p^ = p, Trp = l, p>0, (2) 

i. e. they are Hermitian, have trace one, and are positive semidefinite, i.e. have nonnegative 
eigenvalues. For the sake of simplicity we will often drop the term "semidefinite" in the text 
and write "positive" only, meaning positive semidefinite. All states p satisfy the inequality 
Trp^ < 1, and are classified as pure states (Trp^ = 1) or mixed states (Trp^ < 1). 
The Hilbert-Schmidt distance of two states pi and p2 is given by 

d{pi, P2) := IIpi - P2II := (pi - p2, Pi - 92^^'^ ■ (3) 

Density operators can be represented as a matrix using an orthonormal basis of Ti. The 
standard basis representation is 

pijM ■■= {i3\p\kl) (4) 

where \ij) := \i) ® \j) and {|«j)} is the standard product basis of a (i^ x c^s dimensional 
system. 

A state p is called entangled if it cannot be written as a convex combination of product 
states (l9(], 

p^(j = "^Pip\®p%, (5) 

i 

where p\ and p^ are states of the two subsystems, usually called "Alice" and "Bob" following 
the convention of quantum communication. States a dS]) are called separable, they are not 
entangled and contain classical correlations only. There is still no operational method (i.e., 
no "recipe" ) to decide for a given bipartite state of arbitrary dimension wether it it can be 
written as the convex combination ([5]) or not. This problem is known as the separability 
problem. However, there exist several criteria which help to find the entanglement properties 
of states in special cases, for an overview see, e.g., Refs. |7|, |9|, |20| 



A criterion that is necessary and sufficient, but not operational, is the entanglement 



witness criterion (EWC) 0, S i 

if it can be "witnessed" by some Hermitian operator A, for which 



23l . |2J| . It says that a state p is entangled if and only 



{p,A) = TipA < 0, 
{a, A) =TiaA > 



WaeS, 



(6) 

(7) 



where S denotes the convex and compact set of all separable states, and A is called an 
entanglement witness. We call inequality ([6]) "entanglement condition" and inequality ([7]) 
"separability condition". If there exists a separable state a for which TraA = 0, then A is 
called an optimal entanglement witness. It is closest to the set of separable states and thus 
detects more entangled states than non-optimal witnesses. 

The EWC is a consequence of the Hahn-Banach theorem of functional analysis. It geo- 
metrically corresponds to the fact that an element of a Banach space can always be separated 
by a hyperplane from a convex and compact subset that does not contain the element (see, 
e.g., Ref. and Refs. 0, Q 



for illustrations). Although it is intuitive and simple, the 
EWC is not easy to implement given an arbitrary state p, since in general it is difficult to 
find a suitable witness that satisfies Eq. (jTj), and even more difficult to state that there does 
not exist any witness for this state, which would imply separability. Nevertheless the crite- 
rion plays an important role in the theoretical understanding of entanglement, and has the 
advantage that a witness A corresponds to a physical observable that can be implemented 
in experiments. It therefore allows a detection of entanglement without performing a full 
tomography of the state [2| 
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An operational criterion that is a necessary condition for separability is the positive 
partial transpose (PPT) criterion [Sll- It simply says that a separable state a stays positive 
under partial transposition, 

or^ := (l®r)or > 0, (8) 

where the partial transpose is a transposition T with respect to Bob's system only, pf^i^i : = 
pii^kj- As a proof of Eq. ([8]) one just has to recognize that if we apply the partial transposition 
to a separable state ([5]) the transposition is performed on p^ only, which does not change 
the positivity of the state, and thus the whole separable state stays positive. Therefore, if 
a state p violates the criterion, i.e. it is no longer positive under partial transposition, it 
has to be an entangled state. We call a state that is positive under partial transposition 
PPT, and a state that is not NPT. For dimensions 2x2 and 2 x 3 it can be shown that the 
criterion is necessary and sufficient 
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i.e. any entangled state has to be NPT. In higher 
dimensions, however, there exist entangled states that are PPT, it can be proven that such 
states are bound entangled. Note that the reverse is not necessarily true, it is still an open 
question if all NPT states are free entangled, although there are strong implications that 
NPT bound entangled states exist |32|, |33 . 

In order to find PPT entangled states, one has to employ criteria that are not equal to 
or weaker than the PPT criterion. In principle the EWC is strongest since it detects the 
entanglement of all entangled states, but it is more cumbersome to apply. 

Another useful criterion is the realignment criterion (or cross norm criterion) |3J, |35|, ISy, 



37 



which again is a necessary condition for separability. It states that for any separable 
state the sum of the singular values Si of a realigned density matrix a^ has to be smaller 
than or equal to one, 

~ ^~~ (9) 



E' 



Tr 



^R^R 



< 1 



where {pij,ki)ji '■= Pik.ji- The reahgnment criterion is neither weaker nor stronger than the 
PPT criterion, meaning that it detects some entangled states that the PPT criterion does 
not, and vice versa. Thus an apphcation of both criteria is easy to perform and allows 
a detection of many entangled states, both free and bound entangled, but they do not 
constitute a necessary and sufficient criterion together, since there exist PPT entangled 



states that are not detected by the realignment criterion [37 



III. BLOCK DECOMPOSITIONS AND ENTANGLEMENT WITNESSES 

Bloch decompositions are a convenient way to handle calculations in high dimensional sys- 
tems, since the usage of large matrices can be avoided (for an overview see [38| and references 
therein). There also exist computable separability criteria based on Bloch decompositions 
of states I39I, |40| . 



Let us consider the Hilbert-Schmidt space for a one-particle state first, for example the 
space of Alice's subsystem Aa on TCa of dimension cIa '■= d (all considerations are equivalent 
for Bob). Since the Hilbert-Schmidt space Aa is a vector space of operators, one can 
decompose any element of Aa into a linear combination of operators that form an orthogonal 
basis of the Hilbert-Schmidt space. Let us identify such a basis of d"^ operators with {1, Ai}, 
i = l,...,(i^ — 1. The operators A^ are traceless, Ti Ai = and satisfy the orthogonality 
condition 

TYAiAj = NASij, NaER. (10) 

A one-particle qudit state can then be decomposed into the operator basis as (for example 
a state pA for Alice's subsystem of dimension d) 



PA = -Al + \ \r ' Yn^AA , n, gC, (11) 




where \n\^ = 'Ylii^l^i ^ 1- The coefficient vector n is called Bloch vector, it uniquely 
characterizes the state. The constant ^Jd{d — 1)/Na results from the inequality Tr/?^ < 1. 
The state is pure if and only if |rip = 1. In general not all arbitrary vectors n are Bloch 
vectors, i.e. they do not necessarily imply p > 0, see Remark [H 

A bipartite product state (Jp := pA® Pb on Ti of dimension dA x ds can be written as 
(where Bob's orthogonal basis is {l,_Bj}) 

^ 4-1 4-1 

ap = — (^IdA ® IdB + X] /a ni Ai (g) Idg + ^ fs rrij Id^ ® 5^ + ^ fAfs riimjAi ® Bj 



i=o i,j 



^ <n 1-1^1 1-1^1 f dA{dA-l) , dsidB-l) ,^^. 

nnm,"^/fe e C, |n| < 1, |m| < 1, fA-=\ 77 , IB ■■= \ TT , (12) 



Na ' ' Nf 



where the state is pure if and only if |n| = \m\ 
Any operator O & A can be decomposed as 



4-1 4- 



O = e IdA ® IdB + XI ^' ^' ® IdB + XI ^J ^'iA ^Bj + J2 CijAi (g) Bj , 

4=0 j=0 i,j 

e,ai,bi,Cij e C (13) 



For a given operator O and same dimensions of the subsystems d^ = ds ='■ d, one can 
always find an orthogonal basis in which the coefficient matrix C™^ = (Qj), called correlation 
coefficient matrix, is diagonal: Given an operator decomposition (fT3!) . we have to perform a 
singular value decomposition of C, 

S = UC^W^ , (14) 

where U and V are unitary matrices with entries Uij and Vij and S is the resulting diagonal 
matrix with the d'^ diagonal real positive singular values Sj of C^"^ as diagonal entries. The 
new basis operators Df and Df are then given by a linear combination of the old operators, 

Df = J2<jA„ Df = J2^,,B,, (15) 

j j 

which satisfy the same orthogonality condition, TiDfD^ = N^^ij (and equivalently for Df). 
So we can rewrite Eq. 013p as 

cP-l d^-l 

O = e Id ® Id + J] ri Df ® Id + ^ tj Id (^Df + J2 s^Df ® Df , 

i=0 j=0 i,j 

(16) 

where r^ = ^ • ajUij, tj = ^^ bkV*/^ and Sj = ^ • ^ UijCjiv*i. We call the decomposed operator 
that is written in the optimized way of Eq. fITBl) "singular value optimized" (SVO). Of course 
a product state can then also be decomposed in terms of the new basis, 

d:^-i (f-i 

^ (id ® Id + 5Z /a n, Df ® Id + ^ /ij m, Id ®Df + Y^ JaIb niUijDf ® Df ) , 

1=0 j=0 i,j 

ni,Thj G C , 1^1 < 1 , 1^1 < 1 • (17) 

For our purposes we want to reformulate the separability condition (jTj) of the EWC: 

Corollary 1. An operator C ^ A satisfies TicrC > \/cx E S if and only ifTiapC > for 
all pure product states ap := Pa® Pb- 

Proof. If we have TraC > Vcr G S" then of course also Tro"pC > since the pure pruduct 
states (jp are separable states as well. A separable state ([5]) can be written as a convex 
combination oi pure product states, a = YliiPi^pi because mixed states pA and ps are convex 
combinations of pure states. Thus if Tra^C > 0, it follows that TrcrC = Tr^^ pjCTpC = 
^^PjTrcTpC > since Pi > 0. D 

At first sight Corollary [1] may appear redundant, but it bears the advantage that in order 
to check if a given operator satisfies the separability condition (jTj), we do not have to check 
all separable states but consider pure product states only, which implies a decrease in effort. 
Purity of the states is not essential, but is more convenient in parameterizations. 

For an arbitrary operator basis we use the Bloch decomposition (IT^ to write an Hermitian 
operator C G ^ as 

/ 4-1 4-1 ^ 

C =5\ /ild^ ® IdB + XI ^^ ^^ ® IdB + 5^ ij IdA ®Bj + Y^ CijA, 

\ *=o i=0 i,j 

p := ^/{dA-l)idB-l), 5 G R+ . (18) 



ap 




Note that we are only interested in Hermitian operators C with positive values of 5. If 5 
happens to be negative in the first place we switch to the operator with an additional overall 
minus sign. For (Ia = dB = d we use Eq. (TT6|) to obtain the SVO form 

c=6i {d-i)id(^ Id + Y,^i^f® '^d+Y. ^^ ^'^ ® ^f + 5Z ^'^^ ® ^n ' 

\ i=0 j=o % ) 

(19) 

and get for the expectation value with product states dp from Eqs. ( TT2l) and ( ITSj) (we use cr^ 
in order to conveniently utilize the orthogonality condition flTUl) ) 



which simplifies for dA = ds = dio (using Eqs. ( TT7I) and ( TT9l) ) 



(21) 
Using the above expressions for the expectation values we obtain a condition for TrCfjp > 
in Corollary [T] in terms of Bloch decompositions: 

Corollary 2. Given a decomposition (fT8|) o/an operator C into an arbitrary operator basis, 
the expectation value for any product state (1121) is positive or vanishes, TiCap > 0, if and 
only if 



■' - V^M^^'-"' ^ ^nJt^^''"^^^^''^'"' ^ -' P2) 



J «j 



for all Bloch vectors n, m. For equal dimensions of the subsystems, dA = ds = d, the 
condition ( l22l) can be simplified to 



d s-^ _ ^ / d 



' = V iMi^T) E '>«.• + \/iv:a^3T) Ev*; + E 5.'t'-»: > -i . (23) 

where we used the SVO form fITI?]) of C. 

Proof. The proof is evident from the expressions for the expectation values in Eqs. (120]) and 
(EH). D 



Remark 1. Consider the case when there also exists at least one state p for which TrCp < 0. 
Then C is an entanglement witness if S > —1. Note that by stating "Bloch vector" we mean 
vectors n that correspond to states (i.e. p^ > in Eq. (ITT1) ). For arbitrary dimensions d of 
the Hilbert space an arbitrary vector n G ^C^ ~^ for which p has real eigenvalues does not 
always implicate p > 0, this is only true for d = 2, where the familiar matrix basis out of 
the Pauli matrices or rotations thereof is used. 



Remark 2. It directly follows from Corollaries [T] and [2] that C is an optimal entanglement 
witness if and only if there exists a state p such that TrCp < and Bloch vectors n, m such 
that S = —1. 

Remark 3. For operators C (TT8|) with vanishing coefficients a,, bi, condition (!22|) reduces 
to 

J2^^Jn*m*>-l (24) 

and for cIa = ds condition (I25|) reduces to 

5^5,n*m*>-l. (25) 



This is for example the case if we consider geometric operators (see Sec. HVl) constructed of 
states that are locally maximally mixed, which means their reduced density matrices are the 
maximally mixed states (l/rf^)! and (l/rf^)!. 

Lemma 1. For operators C (TT8|) on a Hilbert space Ti of equal dimensional subsystems, 
dA = ds, with vanishing coefficients di, b^, the expectation value for product states is greater 
or equal to zero, TrcXp > 0, if the singular values Si of the correlation coefficient matrix Cij 
are smaller or equal to one. Si < 1. 

Proof. With vanishing coefficients di, bi, the term S in Eq. fl22|) reduces to 5 = J2ij CijUiUij. 
For dA = ds we can write the operator in SVO form, which gives S = J2i Sin*fh*. With the 
condition Sj < 1 we get 



1^1 



E~ — * — * 
SiU-m- 



< J]s,|n*||m*| < J]|n*||m*| < 1 (26) 



and thus S > —1. D 

Remark 4. Note that Lemma [1] gives only a sufficient condition for satisfying the inequality 
TrcTp > 0. It is necessary for dimensions 2x2 only, since in this case any vectors rTj and fhi 
( TT2I) correspond to states, see Remark [H and with at least one singular value Sj > 1 one can 
easily construct Bloch vectors such that S < —1. For higher dimensions it is possible that 
some Si > 1, and still there exists no Bloch vectors Ui and rhi (that provide p > 0) such that 
S < -1. 



IV. GEOMETRIC ENTANGLEMENT WITNESSES 
Definition 1. A geometric operator G & A is defined as 

G := pi - p2 - {pi, pi - P2)1d , (27) 

where pi and p2 are arbitrary states in A and pi 7^ P2- 



The definition originates from tfie construction of entanglement witnesses in Refs. [23|, |39| , 
with the difference that in our definition the geometric operator (127|) does not yet have to 
be an entanglement witness. The construction (127|) provides TrpiG = and Trp2G' < 0, 

TrpiG = (pi, G) = {pi - Pi, Pi - P2) = , 

Trp2G' = {p2,G) = (p2 -pi,pi -P2) = -||pi -P2f < 0. (28) 

It corresponds to a hyperplane in the Hilbert-Schmidt space A that divides the whole state 
space into states p„ for which Trp„G < and states Pp for which TrppG > 0, see Ref |24l |. 
The hyperplane is orthogonal to pi — p2 since for all states pc on the plane, i.e. that satisfy 
TrpcG = 0, the operator pi— p2 is orthogonal to pc—pi because TrpcG = {pc— Pii Pi~ P2) = 
0. 

Definition 2. A geometric entanglement witness (GEW) Aq is a geometric operator that 
satisfies TrapAc > for all pure product states ap. 



Due to its construction, a geometric entanglement witness (see also Refs. [23|, 124, l39|, 
4ll . l42l . |43| has to witness at least the entanglement of p2. For arbitrary states p2 it is 
easy to construct geometric operators G (Definition [T]) that ensure Trp2G < 0, but difficult 
to confirm that also TrUpG > for all pure product states, which would yield G = Aq- 
Nevertheless, due to their simple geometric construction, geometric operators provide useful 
tools to characterize entanglement, as we will see in the further sections. Other methods to 



construct and optimize entanglement witnesses are given in Refs. [4J, |45|, |46|, |47|, |48 . 

To detect entanglement it is sufficient to consider geometric entanglement witnesses only: 

Lemma 2. Any entangled state is witnessed by a geometric entanglement witness. 

Proof. If p is entangled, then there exists a so-called nearest separable state ctq, i.e. the 
separable state for which the Hilbert-Schmidt distance (^ from p to the set of separable 
states S is minimal, because S is convex and compact. The corresponding geometric operator 
o'o — P ~ {<^o,cro ~ P)^D is an entanglement witness, since the corresponding hyperplane 
includes ctq, is orthogonal to (Tq — p and is therefore tangent to S. For more details on 
nearest separable states see Refs. 23|, |2J, |49| . D 



Geometric entanglement witnesses bear the advantage that they can be "shifted" along 
lines of parameterized states. 

Proposition 1 (Shift method). If a geometric operator 

Gx = Px- p- {px,P\- p)Id (29) 

with a parameterized family of states 

px := Ap + (1 - A)p, 0<A<1, p,pEA (30) 

is an entanglement witness in a parameter region A G [Aj, 1), i.e. if it satisfies TiapGx > 
for all pure product states ap, then px is entangled for A G (Aj, 1]. 

Proof. We consider states px with A, < A < 1 and the geometric entanglement witness 
Ax, = Px,- p- {p\,Px, - p)1d- The expectation value in px is 

TrpA^A, = (pa, ^a,) = {px - px„px, - p) 

= (A, - A)(l - A,)(p - p,p - p) = (A, - A)(l - A,)||p - pf < , (31) 

hence the states pa with Aj < A < 1 are entangled. D 




FIG. 1: Outside- in shift method. On the hne between the entangled state p and the separable state 
a one can detect more entangled states, e.g. bound entangled states, by shifting the geometric 
entanglement witness. 

An effective way to use the shift method of Proposition [T] is to identify p in Eq. fl30p with 
a separable state, and p with a state which is known to be entangled. There are two cases 
where this is of particular interest: 

1. {Outside-in shift.) Starting from the entangled state p, we can detect further entangled 
states along the line in direction to the separable state. By proofing that Gx is an 
entanglement witness for a parameter region Aj < A < 1 (the case A = 1 can be 
included with a suitable normalization of Gx), one can infer that all states px within 
this region are entangled. A reasonable choice for the separable state is the maximally 
mixed state {1/D)1. In this way one can detect bound entangled states, for example 
if we choose a PPT entangled "starting state" p, then we are likely to find more bound 
entangled states along the parameterized line px- The outside-in shift is illustrated in 
Fig. [H See Refs. 42, |43| for application examples. In Ref. [50(| a similar approach with 



parameterized lines between PPT entangled states and the maximally mixed state is 
used to identify families of bound entangled states in the context of robustness of 
entanglement. 

2. [Inside-out shift.) Another application of the shift method is the step-by-step con- 
struction of the convex set of separable states. Here one has to use optimal GEWs 
that correspond to hyperplanes tangent to the set of separable states. Let us assume 
we are given a specific convex subset of states for which we want to determine the 
entanglement properties and that some separable states are known. From these we 
can constuct a kernel polytope of separable states, i.e. the convex hull of the known 
separable states. Then we assign geometric operators to hyperplanes tangent to the 
kernel polytope. For example, an operator corresponding to a plane that includes the 
line between two separable states can be constructed in the following way: Given two 
separable states ai and (T2, the convex line between them is a^ = /xai + (1 — /i)cr2. 
Now we choose an entangled state u, that of course lies outside the kernel polytope, 
such that there exists a /ij with a = a^- for which we have the orthogonality condition 
(o"! — a2,uj — a) = 0. The geometric operator is then given by 

G = a-uj- {a,a-uj)lB (32) 
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and a shift operator Gx between a and u according to Eqs. (129|) and (130|) . The con- 
struction of operators that correspond to boundary planes of the kernel polygon iden- 
tified by more than two separable states is done similarily, using more orthogonality 
conditions and the convex hull between three states. 

Once we assigned geometric operators to the boundary of the kernel polytope, we 
utilize Proposition [1] to "shift" the operators outside and survey the minimum of S 
in Eq. (1221) or (12^ . At one point of the parameterized line (150]) we obtain S = —1; 
the geometric operators become optimal geometric entanglement witnesses. In this 
way we can assemble the shape of the set of separable states for the considered set 
of states and distinguish it from the set of entangled states. It may likely be that we 
have an idea of the shape of the set of separable states that we got from applying 
necessary separability criteria. Then we can use the inside-out shift to verify or falsify 
that shape: The inside-out shifted geometric operators should correspond to optimal 
GEWs when they become tangent to the estimated shape. In this way we get vertices 
of a new polytope, whose boundary planes are shifted again. Thus we either verify the 
estimated shape of separable states, or, if a shifted plane is an optimal GEW before 
it is tangent to the shape, it is an enclosure of all separable states and also entangled 
ones. We require finite steps of this method if the estimated shape is a polygon, and 
(in principle) infinite steps if it is not a polygon, i.e. if it has a curved surface. 

If we have no idea of a possible shape of the set of separable states, or if our estimation 
turned out to be wrong, we can use the inside-out shift to obtain at least a tight 
enclosure polytope. It is a polytope that encloses all separable states but might also 
contain some entangled states, it can be obtained by applying the shift to more than 
one kernel polytope. Both situations are sketched in Fig. [2J 

The difficult part of Proposition [1] is to prove that G\ is an entanglement witness, in partic- 
ular the verification of the separability condition ([7]) . To accomplish this we can efficiently 
use the previous corollaries and lemmas, which will be demonstrated by the example of the 
next section. 



V. ENTANGLEMENT PROPERTIES OF A FAMILY OF TWO-QUTRIT 

STATES 

An interesting set of states is the magic simplex of two-qudit states (dimension d x d) 
5ll . |52|, |53|. It is the set of all states that are mixtures of Bell states Pnm, 



W := < ^ qnmPnm | ^nm > 0, ^ g„m = 1 |- , (33) 

[^n,m=Q n,m ) 

where the d^ operators (the Bell states) 

Pnm := (f/„^®l)|0+)(0+|(t/L®l) (34) 

form an orthogonal basis of the d x d dimensional Hilbert space and the vector state |0^) 
denotes the maximally entangled state 

k^) = ^Y.\3)^\3)- (35) 

j=0 
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FIG. 2: Left: Sketch of the inside-out shift with an estimate of the shape of separable states, which 
in this case coincides with the true set of separable states, pictured by the circle. We start with a 
kernel polytope (black triangle) and shift the boundary planes outside until they become optimal 
GEWs, which are tangents to the circle (dashed lines). In this way we can draw a new polytope 
(hexagon, grey). In the next steps (not illustrated) the boundaries of the new polytope are shifted 
and we gain a new polytope, and so on. In this way we reconstruct the circle shape. 
Right: Sketch of the inside-out shift where we do not rely on an estimate of the shape of separable 
states. The true set of separable states is again pictured by a circle. Here we get a first enclosure 
polytope (biggest triangle with dashed lines), by shifting the boundaries of a first kernel polytope 
outside (dark grey triangle). A tighter enclosure polytope (hexagon with dashed lines) is obtained 
by shifting the boundaries of a second kernel polytope (small black triangle) outside. The light 
grey areas mark states inside the enclosure polytope that are not separable and thus account for 
the deviation of the enclosure polytope from the true set of separable states. 

The unitary operates Unm are the Weyl operators 



U„ 



d-l 

E 

fc=0 



e d 



kn 



\k){{k + m) mod(i| 



(36) 



which have been introduced in classical theories of discrete phase space and appear in quan- 



tum information theory in the context of quantum teleportation [5J]. In the teleportation 
protocol the Bell state basis flMl) is the higher dimensional generalization of the two-qubit 
basis and the Weyl operators Unm are the analogue of the Pauli operators, they correspond 
to the operators Bob has to apply in order to obtain the teleported state. The reduced 
density operators of states that are elements of the magic simplex are maximally mixed, but 
not all two-qudit states with maximally mixed reduced density operators are elements of the 
magic simplex, apart from dimension 2x2, where all locally maximally mixed states are 
included in the tetrahedron of all Bell state mixtures [53] . Furthermore the magic simplex 
has a high symmetry in the phase space of the coefficients nm, for a detailed discussion see 
Refs. Sir^ 



1 sy i 

Hi. 



The Weyl operators fl36l) form an orthogonal operator basis. 



Trf/Lf// 



ij 



d 5ni 5; 



mj 



(37) 



and hence can be used for Bloch decompositions. The Bell states Pnm (13^ can be expressed 
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with Weyl operators as (where the indices have to be taken mod d) 

d-l d-1 

^^^ = ^11 e'5'('='-"'") Ui^ ® U^i^ = ^ 5Z c,^ Ui^ ® U^ln^ ■ (38) 

m,l=0 l,Tn=0 

Obviously the Bloch vectors corresponding to the Bell states have a diagonal but in general 
complex coefficient matrix (qj), where i counts the different combinations of /m and j those 
of kn, and Cu := q^- The singular values of the coefficient matrix are Sj = |cij| = \cim\- 

Note that a construction of the type (1M|) can be done with any unitary operators that form 
a matrix basis of the Hilbert-Schmidt space, obtaining other bases of orthogonal maximally 
entangled states. 

A subset of the magic simplex of two-qutrit states (dimension 3x3) that reveals interesting 
entanglement characteristics is the three-parameter family |42l . |43| 

1 — <y — 3 — ^ d 1 

Pa,/3,7 := ^ ^1 + «^oo + I {Pm + P20) + ^ (Poi + Pii + P21) , (39) 

where the parameters are constrained by the positivity requirement Pa.p,-^ > 0, 

7 
«< 2/^ + 1-7, a<-/3 + l-7, 

«<-/3 + l + 27, a>^-i + i7- (40) 

The family of states (!39|l contains a one-parameter family of states that have three entangle- 
ment properties; they can be separable, PPT entangled and NPT entangled. We call them 
Horodecki states ph J16l] . 

P6=^|0'+>(0^| + ^a+ + ^(7_, 0<6<5, (41) 

and, according to our parametrization, 

6-6 ^ 2h 5-26 

Pb ■■= Pa,f3,^ with a = -^, (3 = - — , 7 = — — . (42) 

Using the PPT criterion we find regions of PPT and NPT Horodecki states: They are NPT 
for < 6 < 1, PPT for 1 < 6 < 4 and again NPT for 4 < 6 < 5 . In Ref. ^ it is shown 
that the states are separable for 2 < 6 < 3 and bound entangled for 3 < 6 < 4. 

Now let us apply the PPT criterion ([S]) and the realignment criterion i^ to our three- 
parameter family fl39l) . The PPT criterion provides the following parameter constraints for 

PPT states Pa,l3,'y'- 

r. 1 1 

a < -3 h -7 , 

- ^ 2 2 



a 



< l(^-2 + ll/3 + 3v^) , a>^(-2 + ll/3-3v^) , (43) 



where A = 4 + 9/5^ -|- 47 — 77^ — 6/5(2 + 7). Hence all states pa,p,t with constraints fj43l) are 
either bound entangled our separable, whereas the others are NPT entangled. 
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From the realignment criterion we obtain the constraints 



«< 4(6 + 11/5-7 -Ai) (44) 

« < 4 (6 + 11/? - 7 + Ai) (45) 
lo 

«> ^(-6 + 11/3-7 -A2) (46) 

« > 4 (-6 + 11/5 - 7 + A2) (47) 
lo 



where 



Ai := ^4 + 36/? + 81/?2 - 127 - 54/^7 + 8872 and 

A2 := ^4 - 86/5 + 8I/32 + 127 - 54/^7 + 8872. (48) 

Only constraint flH|) is violated by some PPT states, which thus have to be bound entangled. 
The PPT entangled states exposed by the realignment criterion are therefore concentrated 
in the region confined by the constraints 

«<^/5 + l-7, a<^(-2 + ll/5 + 8v^), a> ^ (6 + 11/3 - 7 - Ai) . (49) 



The three-parameter family (1391) also bears the advantage that it can be nicely illustrated by 
the Euclidean geometry. To do this, note that the orthogonality conditions of the Hilbert- 
Schmidt space A have to be transferred correctly, which is achieved by choosing a nonorthog- 
onal and differently scaled coordinate system of parameter axes a, (3, and 7. They are chosen 
such that they each become orthogonal to one of the boundary planes of the set of the three- 
parameter family of states, given by the positivity constraints (HOj) . In order to calculate 
quantities well known in an Euclidean space spanned by an orthogonal equally scaled coor- 
dinate system, we have to transform points of the non-orthogonal coordinates (a,/3,7) into 
points of orthogonal coordinates (a, b, c) and vice versa by 

« = «-^/5-^7, 6 = ^(8/3-7), c=^7. (50) 

In Fig. [3] the three-parameter family of states Pa,/3,7 ( !39l) including NPT entangled, PPT en- 
tangled (bound entaiigled) and further PPT states are illustrated in the Euchdean geometry 



picture. In Refs. [38|, |42j we applied the outside-in shift method to detect most of the bound 
entangled states (1491) . where a version of Lemma [T] (for Weyl operator decompositions) was 
used to show that for particular parameter regions geometric operators correspond to geo- 
metric entanglement witnesses. The geometric shifting operators Gx fl29l) were constructed 
on lines between bound entangled starting states p = pf^ on the Horodecki line (j^Tl) and 
the maximally mixed state, p = (1/9)1. 

Actually all PPT entangled states of Eq. fH9l) . Fig. [3l can be detected using Lemma [H 
To see this, we construct tangent planes onto the surface of the function 

a= ^ (6 + 11/3-7 -Ai) (51) 
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FIG. 3: Illustration in of the family of states Pa,p,-y (|39|) in the Euclidean geometry. Left: All 
states Pa,/3,'r lie within a pyramid due to the positivity constraints (|10|) . The dot represents the 
origin of the coordinate axes, which is the maximally mixed state (1/9)1. The line (blue) on 
the left boundary plane represents the Horodecki states (1411) . Right: Illustration of the PPT and 
realignment criteria. The cone with tip on the right vertex line of the pyramid contains the PPT 
states, which is intersected by a cone (tip on the left boundary plane of the pyramid) of states that 
satisfy the realignment criterion, and hence PPT entangled states can be revealed (translucent 
yellow region). All other states of the pyramid are NPT entangled. 



from the realignment criterion (jHj), where we use orthogonal coordinates (!50|) . In this way 
we can assign geometric operators to the tangential plane by choosing points a inside the 
planes and points b outside the planes such that a — b is orthogonal to the planes. Since 
the Euclidean geometry of our picture is isomorphic to the Hilbert-Schmidt geometry, the 
points a and b correspond to states Pa and pb and we can construct the geometric operator 
accordingly, 

Grc = Pa- Pb- (pa, Pa " Pb) I9 • (52) 

This operators ( l52l) are linear combinations of the three-parameter states Pa,i3,j which are 
linear combinations of the Bell states Pnm (1341) and can be written as a Bloch decomposition 
using Eq. (!38l) . First we need to define some expressions of Weyl operator combinations. 



Ui := f/oi ® f/oi + Uo2 Uo2 + Uu ® t/-ii + f/12 ® t^-12 + U21 f/-2i + U: 



22 



f/_ 



22 



f/2 := u' + a 



II 



with f// := U- 



10 



f/_ 



10 



a 



II 



u- 



20 



f/_ 



20 



(53) 



The geometric operators fl52l) corresponding to tangent planes in points (at,/3t,7t), where 
at is a function of (3t and 7t, given by the realignment function flSTj) . are 



Gr 



a{21-Ui + cUi + c*U. 



*ttII\ 



with 



- (-2 - 9A + 37t + 3A,) , 
6b 

97t^ + (-2 - 9A + 37t)Ac + V37t (2 + 9A - 37t + 3A,) 

(2 + 9A)2 - 6(2 + 9A)7t + 367t' 

:= V4 + 36 + 81A' - 127t - 54A7t + 337? ■ 



(54) 



15 




® ® 




FIG. 4: Left: The entanglement properties of the three-parameter family on the boundary plane 
(j56p where the Horodecki states are located. The triangular region (green) contains the separable 
states, the bound entangled states are located in the parabolic region (yellow), and the remaining 
states are NPT entangled. Right: The kernel polytope is a polygon (green) that includes states 
that are necessarily separable. 



The singular values of the correlation coefficient matrix are the absolute values of the coef- 
ficients —1, c and c* in Eq. fl5^ . which are all one. 



{s,} = {1,1, 1,1, 1,1, 1,1} 



(55) 



and therefore, according to Lemma [H the geometric operators G^e are entanglement wit- 
nesses that detect the entanglement of all states "above" the corresponding planes, thus also 
the bound entangled states in the region of Eq. (B^ . 

We might ask ourselves if the PPT entanglement of Eq. fH^ . revealed by the realignment 
criterion and also by GEWs, is all there is for the three-parameter states ([31]). Or, to put it 
differently, are all the three-parameter states that satisfy both the PPT and the realignment 
criterion separable? We can answer this question by using GEWs and the inside-out shift 
method. The entanglement properties of the states on the boundary plane 



7 



a 



/3 + 1-7 



(56) 



of the positivity pyramid are already fixed. The realignment function flSTj) and also the 
GEWs Gre dSH) draw a triangle on this plane, whose vertices are separable states. The 
tip of the triangle is a separable state since it is PPT and_7 = (all PPT states of the 
two-parameter subset 7 = are separable, shown in Ref. [5l|), the other two, at 7 = 1 
and ' 1 = —1 are simple mixtures of Bell states Pnm that are also shown to be separable in 
Ref. 5l| • So the the triangle is the convex hull of the three separable states and thus has 
to be separable. For an illustration of the entanglement properties on the boundary plane 
dMD see Fig. H 

But what about all the three-parameter states (!39l) ? First, we construct a kernel polytope 
of those states that are necessarily separable. This can be done by identifying five separable 
states that serve as vertices for the kernel polytope. Three arise from the two-parameter 
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subset 7 = 0, where all PPT states are separable, the remaining two vertices are the separable 
states with 7 = 1 and 7 = — 1 on the boundary plane (1561) . The resulting kernel polytope 
is a polygon with five vertices, see Fig. |H Alternatively, one can also use sufficient criteria 
for separability to construct a kernel polytope of separable states. In Ref. [55| a sufficient 
separability criterion is presented that is shown to be applicable for states of the magic 
simplex. 

In Ref. J43|] it remained open if this polygon contains all separable states of the three- 
parameter family fl39|) . which would imply much greater regions of bound entanglement than 
detected before. Here we want to show that this is not the case. 

We can assign geometric operators to four boundary planes of the kernel polygon, in 
the same way as we did for the planes on the realignment surface, see Eq. f l52|) . where we 
use the geometric isomorphism again. We call the four geometric operators G'!^,G'^, which 
correspond to the following planes given by three vertex points (see Fig. H]): G" to 1,3,4, 
G" to 3,4,5, Gi^ to 1,2,3, and G'^_ to 2,3,5. The plus and minus sign indicates the side 
with positive or negative values of the parameter 7. The operators are 

Gl = a{21-Ui + cUi + c*Ul^) , with a = — , c = -1 ± VSi 

63 

G'i = a {21 + Ui + cUi + c*Ui^), with a = ^, c = -1 ± VSz (57) 

The boundary planes can be easily shifted along parameterized lines through their normal 
vectors, and so can the assigned geometric operators (1571) . Note that we have a simplified 
picture of locally maximally mixed states, see Remark [3l The operators G!^, G^ themselves 
are not entanglement witnesses, since the condition (IMl) can be numerically shown to be 
violated (see Corollary [2]). The singular values are again the absolute values of the correlation 
coefficients, {sj} = {1, 1, 1, 1, 1, 1, 2, 2}, hence Lemma[T]does not give an answer. It is difficult 
to show a violation analytically because of the complex Bloch vector geometry of qutrits, 
see Remarks [1] and HI In order to check the condition flM|) we minimize the left-hand term 
numerically by varying the possible Bloch vectors n*,m*, restricted by the condition p > 
with p of Eq. ( ITTl) . Shifting the operators outside, we find a minimum S" = — 1 of condition 
([21]) when the planes become tangent to the shape enclosed by the PPT and realignment 
criterion, achieving new vertices at the touch points. Employing the inside-out shift method, 
see Sec. |IV] and Fig. |2l we construct a new polygon with the new vertices, and assign new 
geometric operators corresponding to the new boundary planes. Shifting the new operators 
outside, we again find the minimum S* = — 1 at planes tangent to the PPT and realignment 
shape. Therefore there is a very strong implication that the PPT and realignment shape, 
seen as the two-cone shape in Fig. [3l is the shape of the separable states. Fig. ([3]) thus is a 
picture of all entanglement properties of the three-parameter family. 



VI. SUMMARY AND CONCLUSION 

We use the concept of Bloch decompositions and entanglement witnesses to detect the 
entanglement properties of arbitrary dimensional bipartite quantum states. In particular 
we show how to reformulate the conditions of the entanglement witness criterion by using 
Bloch decompositions (Corollary [2]) and formulated a sufficient condition for an operator to 
be an entanglement witness (Lemma [T]). 
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We give the definition of a geometric operator and a geometric entanglement witness 
and explain two methods of "shifting" it (Proposition [1]): One for the detection of bound 
entangled states, the outside-in shift, and one for the detection of separable states and for 
the construction of the shape of the set of separable states, the inside-out shift. 

Finally we apply the previous results on a family of three-parameter two-qutrit states 
that are part of a simplex in the state space of two qutrits, the magic simplex. We show 
how to detect bound entangled states and construct the shape of separable states for this 
family. The results can be conveniently illustrated by the Euclidean geometry. 

Our approach to entanglement detection is guided by the geometrically intuitive way 
of using entanglement witnesses. The construction of geometric entanglement witnesses 
directly uses the fact that entanglement witnesses correspond to hyperplanes in the Hilbert- 
Schmidt geometry. In this way it becomes easier to apply geometric operations like the 
shifting of planes. Using Bloch vector decompositions of operators and states we can fur- 
thermore simplify the conditions that have to be satisfied such that a geometric operator 
is a geometric entanglement witnesses. The construction of geometric entanglement wit- 
nesses does not rely on special properties of the states, i.e. it can be done for NPT or PPT 
entangled states likewise. 

The presented example is relevant in many aspects. First of all the states of the magic 
simplex are a higher dimensional analogy of Bell-state mixtures of the two-qubit case that 
are relevant for quantum communication tasks, as explained in Sec. |Vl Furthermore it is 
interesting and surprising that this particular three-parameter family includes the Horodecki 
states that were among the first examples of bound entangled states. Thus the three- 
parameter family can be viewed as a more-parameter extension of the Horodecki states that 
includes even more bound entangled states. Finally the three-parameter states allow a nice 
Euclidean illustration that makes the regions of entangled, bound entangled and separable 
states visible. 

Throughout the paper we restrict ourselves to bipartite states. Of course a multipartite 
extension is trivially possible if we only want to distinguish between states that contain 
entangled states in any of its particles and states that are fully separable into all particles. 
The definition of separable states just has to be extended with additional tensor products 
respectively. In the case of multipartite states one can distinguish between the distillability 
of states into entangled states of a fixed number of particles |56|. Entangled multipartite 
states can themselves be classified in different ways, for example with respect to the number 
of particles that are entangled. For details see, e.g., Refs. J57l. ISSl. 1591 . 
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